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I Abstract. Based on random matrix theory, we compute the hkehhood of saddles and minima 

^^ in a class of random potentials that are softly bounded from above and below, as required 

^ for the validity of low energy effective theories. Imposing this bound leads to a random 

mass matrix with non-zero mean of its entries. If the dimensionality of field-space is large, 
inflation is rare, taking place near a saddle point (if at all), since saddles are more likely 
^r^' than minima or maxima for common values of the potential. Due to the boundedness of the 

-Y-\ potential, the latter become more ubiquitous for rare low/large values respectively. Based on 

(O the observation of a positive cosmological constant, we conclude that the dimensionality of 

^^ field-space after (and most likely during) inflation has to be low if no anthropic arguments 

. . are invoked, since the alternative, encountering a metastable deSitter vacuum by chance, is 

. 5^ extremely unlikely. 
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1 Introduction 

What is the best model of inflation? Whether or not a model is natural, compelling or 
beautiful is often in the eye of the beholder, but ultimately consistency with observations is 
the deciding factor. If a single observation is at odds with a model, it should be abandoned. 

Models of inflation in string theory are varied, see [1-6] for reviews, involving moduli 
fields, branes and fluxes among other ingredients. One of the first well studied models in 
string theory is the KKLMMT [7] brane inflation [8-15] setup, where all but one effective 
degree of freedom are carefully stabilized. However, soon it was realized that in the absence of 
a compelling a priori reason, several dynamical degrees of freedom appear naturally, leading 
to a plethora of multi-field inflationary models, see i.e. [16-30] for a small, non-representative 
selection. 

Unfortunately, predictions for correlation functions in the cosmic-microwave-background 
radiation are often degenerate, offering little discriminating power [31-37] so far. Non- 
Gaussianities (NG) can provide such a discriminator if it were observed, since simple single- 
field models are unable to generate them [38]; however, there is no compelling reason why 
NG should be large in multi-field models either [39, 40]: to (temporarily) boost non-linearity 
parameters, sudden events are required, such as a sharp turn of the trajectory, a sudden 
(temporary) increase of slow roll parameters, a decrease of the speed of sound, or particle 
production events, to name a few^. Although such events can occur, and they might be 
natural for certain model classes, it is not clear how generic they are. 



^To make accurate predictions for i.e. fluctuations in tlie CMBR, one has to follow perturbations until an 
adiabatic (i.e. single-field) regime is reached [41, 42]. 



This uncertainty is compounded by our limited understanding of the landscape in string 
theory [43]. We know little about the origin of the low-energy, effective-field-theory potential 
in our universe; even if we understood all of the relevant physics, the potential would most 
likely be far from simple, involving many parameters whose values we would be unable to 
compute from first principles. An example is the Denev-Douglas landscape [44], which was 
modelled recently as a random potential [45] (see also [46]). 

Fortunately, if a potential of dimensionality -D ^ 1 is sufficiently complex, some prop- 
erties, such as the distribution of minima or saddle points, can be predicted quite accurately 
using random matrix theory (see i.e. [47, 48] for a textbook introduction) based on the feature 
of universality [49-53] . 

We follow this approach with the aim to answer a simple question: 

What is the inflatons ' most likely resting place, if inflation is driven by multiple fields in a 
random potential that is softly bounded from above and below, with equal likelihood for 

positive and negative values? 

The quantitative answer needs to be compared with the observation of a small, but positive 
cosmological constant (CC) in our universe. 

The novel feature in our analysis is the incorporation of upper and lower bounds, as 
needed for the viability of effective-field-theory, by means of a Gaussian suppression. As 
the dimensionality of field-space increases, we find an extraordinarily strong preference of 
negative values; consequently, almost all inflationary trajectories end up in an anti-deSitter 
(ADS) minimum, at odds with observations. In this sense, multi-field infiationary models of 
this type are strongly disfavored, the more so the higher the dimensionality of field-space D 
{D ~ 500 on the string landscape). While completing this study, a paper by Chen.et.al. [46] 
came to similar conclusions, but the full ramifications for multi-field models were not spelled 
out. 

One could resort to anthropic reasoning to reconcile this inconvenient observation with 
multi-field models - after all, eternal infiation appears to be able to populate all vacua in a 
landscape [54] and one of them should be just right for us^. We refrain from such a reasoning 
throughout most of this paper, primarily due to our personal bias: anthropic reasoning 
had been invoked repeatedly in science in the past, oftentimes precluding much simpler 
explanations. Further, where it is applicable, anthropic reasoning seldom explains anything 
beyond the observation it was invoked for in the first place - it primarily puts one's mind 
at ease. Weinberg's prediction of an upper bound on a positive CC based on anthropics 
(merely requiring the formation of gravitationally bound objects) consistent with current 
observations is a notable and admittedly relevant exception [55]. Similar lower bounds for 
a negative CC exist [56], which are, on the one hand, weaker since a negative CC actually 
enhances the build-up of structures (on the other hand, a negative CC leads to a re-collapse 
of the universe) . 

With this caveat in mind, we take as an agnostic working hypothesis that the CC is not 
entirely set by anthropics. If it turns out that almost all conceivable models are inconsistent 
with a positive CC, one may revert back to anthropic reasoning, see Sec. 5.2, but it appears 
to us that there is no lack of (single or few-field) infiationary setups with the ability to 



^Note that this hne of reasoning opens up a whole can of worms with regard to the measure, which is ill 
defined in an eternally inflating Universe [57]. 



accommodate deSitter vacua^ that avoid the additional challenge brought forth by a large D 
during inflation. 

In addition to drawing conclusions about the final state after inflation (theoretically 
and numerically) , we investigate the feasibility of inflation in random potentials in the first 
place, following prior research by Tegmark [58] as well as Frazer and Liddle [59, 60]. We 
find that inflation is unlikely (the more so the higher D is), and if it occurs, it almost always 
takes place near a saddle point. Further, we observe that a short burst of inflation is more 
likely than a longer one, P oc N^^, where A^ is the number of efolds and /3 a model specific 
parameter of order /3 ~ 3, based on the numerical runs in a particularly simple (from a 
computational point of view) random potential [59, 60], that retains the feature of being 
softly bound from above and below, recovering and extending results by Frazer and Liddle 
[60]; this result is consistent with the theoretical prediction (3 = 3 for inflation at a saddle 
point that was reached in [62, 63]. In essence, we provide additional theoretical explanations 
for some of the observations made in [60]. We would like to point out that many of the 
theoretical conclusions were reached before, albeit scattered over several publications and 
years ^. 

The outline of this paper is as follows: to investigate the distribution of minima, maxima 
and saddle points, we first provide a brief review of known results in random matrix theory. 
Sec. 2. Random matrix theory is needed, because the Hessian is, to a good approximation, a 
random matrix if enough random parameters enter the low energy effective potential. As a 
concrete test-bed, we employ a simple random potential retaining the feature of boundedness 
in Sec. 3. This potential is chosen to be relatively smooth to support inflation. Aided by 
this toy model, we draw general conclusions about the final state after inflation. Sec. 3.1 and 
Sec. 3. 2, based on universality in random matrix theory. We compare our analytic estimates 
with numerics in Sec. 3.2, paying particular attention to the dependence of the minimas 
probability distribution on the dimensionality of field-space. We investigate the feasibility 
of inflation in Sec. 4, again comparing theoretical expectations with numerical results in the 
toy-model. The consequences for the final state after inflation and the resulting unlikeliness 
of multi-field inflation are elaborated on in Sec. 5; in Sec. 5.2 we comment on a possible 
way to avoid these conclusions if anthropic arguments in an eternally inflating universe are 
invoked. We conclude in Sec. 6. 

2 Some Properties of Random Potentials and Inflation therein 

Here, we present some known results based on random matrix theory in the realm of infla- 
tionary cosmology. Readers familiar with this topic may jump to Sec. 3. 

2.1 Likelihood of Saddle- Points VS. Maxima/Minima 

In [45, 61] random potentials with multiple degrees of freedom and their inflationary tra- 
jectories were studied with the aim to extract statistical properties of the distribution of 
minima, maxima and saddle-points (see also [46]), as well as the feasibility of inflation on 

^We do not aim to address the cosmological constant problem here, that is why the magnitude of the CC 
is so small compared to predictions in the standard model of particle physics; we merely point out that the 
prediction in a class of multi-field inflationary models is a large negative contribution to the CC. 

** Random matrix theory is a well studied branch of mathematics with many applications in i.e. solid state 
physics; one of the first applications in cosmology was in the less known, but excellent paper by Azami and 
Esther [61], which has recently caught the attention of Marsh et.al. [45] and Chen et.al. [46]. The preference 
of shorter bursts of inflation near a saddle point over longer ones goes back to [62] (see also [63]). 



such a landscape. Here, we would like to briefly review some of their main results, before 
extending them. For a general introduction to random matrix theory see i.e. [47, 48]. 

Consider D scalar fields (/?j, z = !,...,£•, with potential V{(pi, . . . , ifn)^ canonical kinetic 
terms and a flat fleld-space metric as a model for the landscape in string theory [43]. If D 
is large, a box with volume (2A) with A < Mp is expected to possess a large number of 
critical points with W = (here V is the gradient in fleld-space), say a^ with a ~ few, and 
a smaller but still large number of extrema^. The proportion of extrema to saddle points 
is important for several reasons: flrst, anthropic arguments to explain the smallness of the 
cosmological constant rely on the presence of many minima (~ 10^*^^) besides a dynamical 
mechanism to scan this landscape. Second, to drive small fleld inflation, a close encounter 
with a saddle-point is needed before the fields get trapped in a minimum with small but 
positive energy (large field models suffer from the r] problem, that is quantum corrections 
usually spoil the needed fiatness for super-Planckian field excursions (such corrections are 
one of the main motivations to consider random potentials); further, inflation driven in a 
metastable minimum leads to an observationally ruled out open universe after tunnelling in 
the absence of a subsequent slow roll phase. Thus, one should expect that inflation near a 
saddle point is also the most likely occurrence as D is increased in a random potential.) 

Extrema are classifled by the Hessian, or mass-matrix. 



H 









(2.1) 



and if y is a random potential, the Hessian becomes a random matrix. However, in general 
the entries are not devoid of structure since the potential satisfies constraints; for example, in 
the super-gravity models of Denef and Douglas [44] , both the Kaehler- and superpotential are 
"random" , leading to a Hessian that can be approximated as a sum of Wigner and Wishart 
matrices [45] (see also [46]). In [61], a simpler case was studied from a pure phenomenological 
viewpoint, assuming the Hessian to be a real Wigner matrix (symmetric, Gaussian distributed 
entries with identical variance and zero mean). Given a Hessian, one needs to find the 
probability for H to be positive definite: if the Hessian has only positive/negative eigenvalues, 
a local minimum/maximum is present, and for mixed eigenvalues the critical point is a saddle 
point. 

To compute the probability of a minimum, one needs to derive the joint probability 
distribution of eigenvalues, followed by an investigation of fiuctuations of eigenvalues to 
compute the probability that the lowest eigenvalue is positive. Performing this task for small 
fiuctuations of the smallest /largest eigenvalues goes back to Tracy and Widom [64] (see also 
[65], and [66, 67] for extensions), while large fiuctuations, which are needed in our case, 
have been the subject of a series of papers by Majumdar et.al. [68-73], see [45] for a brief 
review. The eigenvalue spectrum and fluctuations of extreme eigenvalues turn out to become 
independent of the statistical properties of the matrix elements in the large D limit [49-52] , 
and it appears that this universality still holds if entries are highly correlated [53] . Motivated 
by this property, we focus on simple toy models in this paper, anticipating that conclusions 
are true more generally. 



^By considering a finite box, we are implicitly imposing a measure on the landscape; computing the 
frequency of saddles to minima/maxima is ill defined on an infinite field-space without imposing further 
restrictions. 



Let us start with the setup of [61], focussing on two extreme cases: first, assume that 
cross couplings between the fields are small so that the Hessian at the critical point can be 
diagonalized by a field-redefinition. Second, assume that cross-couplings are large so that 
the Hessian is a random matrix with elements of comparable magnitude. 

Near a critical point the potential can be written as 

i i^j 

where the fi each have at least one extremum {dfi/d(pi = 0), and eijipiipj parametrize cross- 
couplings. In the first (unlikely) case where the cross couplings, Cij, are small compared to 
the diagonal elements in H, the Hessian is diagonal to a high degree. The eigenvalues are 
then independent random variables with zero mean and some variance. Due to universality, 
the matrix elements are taken to be independent Gaussian variables in [61]; the probability 
of a minimum/maximum becomes 

P=^- (2.3) 

In other words, the likelihood of a large fluctuation shifting all the eigenvalues to the same 
side of zero is 1/2 . 

If, on the other hand, the cross couplings are of the same order as the diagonal elements, 
the joint probability distribution obeys Wigners semi-circle law for the Gaussian ensemble 
[47] and the probability of a minimum/maximum scales as 

P oc e-^^" , (2.4) 

with constants c ~ 1/4 and p = 2 [61]. Intermediate e^j interpolate between the two extremes; 
for instance, decomposing H = A + B, where A is diagonal with variance a a of its Gaussian 
distributed entries and B random with variance as, it was found numerically that [46] 

P oc e-^^^'-''^ (2.5) 

where 

C2 ^ 0.000395y + l.OSy^ - 2.39y^ + ©(/) (2.6) 

C2/C1 « 0.0120 + 2.99y - 12.2y2 + 1650y^ + ©(/) (2.7) 

and y = aA/c^B "^ 0.1. In general, the constants also depend on the joint probability 
distribution of the eigenvalues. For complex Wigner matrices one finds c = ln{3)/2 and 
p = 2 [61, 68, 70], while the more complicated Denev-Douglas landscape leads to c ~ 0.08 
and p ~ 1.3 [45]. 

The main lesson is that large fluctuations of eigenvalues, as needed for a positive definite 
Hessian, are extremely rare if the joint probability distribution is reminiscent of Wigner's 
semi-circle law, much more so than the naive intuition based on statistically independent 
eigenvalues suggests, exp(— cD^) with 1 < p < 2 vs. 1/2^. Consequently, almost all critical 
points are saddle points. 

An intuitive visualization of the joint probability density of the Wigner ensemble in 
terms of a Coulomb gas is due to Dyson [74, 75]: the eigenvalues correspond to the positions of 
D charged beads undergoing Brownian motion on an infinite string, with a confining quadratic 
potential at the origin and logarithmic repulsion between the beads. The probability of 



finding only positive/negative eigenvalues corresponds to the probability of finding all beads 
to the left /right of the origin. To shift the outermost bead, a displacement of order 0{yD) 
is needed. Since the potential is quadratic, such a displacement comes at an energy cost of 
order 0{D); the large fluctuation above requires ~ D/2 beads to shift, so that the statistical 
cost is of order ~ exp(— D^) [45], in agreement with [61, 68, 70]. 

Applied to inflationary cosmology, we can deduce that inflation is most likely driven 
by fields near a saddle point, not by fields trapped in a metastable vacuum, or even more 
unlikely, via chain inflation models [76, 77] (see [78] for a recent (critical) feasibility study). 
This expectation was confirmed in the KKLMMT setup investigated in [63], as well as for 
the random potentials of [59, 60, 79]. 

Another important application is the search for deSitter vacua in string theory that 
are compatible with a small cosmological constant: uplifted AdS- vacua in type II A string 
theory acquire off-diagonal terms in the mass matrix, often pushing the lowest eigenvalue to 
negative values and thus destabilizing the vacuum [46] . 

However, one crucial feature is missing in the above arguments: effective potentials 
used for inflation, as the ones in [59, 60, 63], are only valid at low energy \V\ <C Mp] thus, 
by construction, the potentials have to be bounded from above and below to remain in the 
validity range of effective-field-theory; as a consequence, minima are more likely close to the 
lower boundary and maxima close to the upper one. The inclusion of this feature is the main 
goal of this paper, see Sec. 3; we comment on consequences for inflation in Sec. 4. 

2.2 Likelihood of connecting two Saddle-Points during Inflation 

If two or more periods of inflation are to be connected during the last sixty efolds of inflation, 
as envisioned in folded inflation [80], it is necessary that the end of an inflationary trajectory 
connects to another saddle-point within a distance H inside field-space. The number of 
extrema present in the cubic toy-landscape (2A)^ with A ~ Mp introduced in Sec. 2.1 falls 
below one if a < 2Z)^' ^, so that the separation r between adjacent minima is too large. For 
larger a the likelihood of finding two adjacent inflationary saddle-points can be estimated as 
[61] 
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under the assumption that the energy scale on field-space is reduced compared to the Planck 
scale by l/\/D, H ~ M"^ /Mp ~ 1/ D. This probability decreases sharply as the dimensional- 
ity of field space increases. Thus, Azami and Easther concluded that inflation is more likely 
to happen near a single saddle-point, if at all. Again, this is in line with the recent case 
studies in [59, 60, 63]. 

2.3 Likelihood of obtaining N e-folds of Inflation near a Saddle-Point 

We saw in the last section that inflation on a random potential most likely occurs at or near 
a single saddle-point. Since, by construction, saddle-points lie within shallow regions, the 
inflationary trajectory can often be approximated by a straight line in field space (an effective 
single-field model can be used), traversing an inflection point. 

Let us therefore follow [62, 63], where inflation driven by a single field occurs near an 
inflection point, located at a position where V" = and V is small. The origin of y? is chosen 
to correspond to the zero of V'\ such that 

V{^)^co + ciip + c^^^ + ..., (2.9) 



and the Cj are constants. In the regime where the constant term dominates and the q are 
small, the number of e-folds of inflation is of order 

iVe~-^xO(l). (2.10) 

To obtain the probability of N^. e-folds of inflation, Freivogel et.al. [62] suggest to compute 

P(iVe) = j IlUdiiF{ii,...i.)5{N,- f{i^,...i,)) , (2.11) 

where ^j are the model's parameters, / is the number of e-folds as a function of ^j, and F 
is a measure of the parameter space. Taking a flat measure, F = (cq, ci, C3) = 1, Eq. (2.11) 
simplifies to 

P{Ne) = j dcidc^5{Ne - -==) . (2.12) 

To evaluate the integral, one uses / f{x)5{g{x)) = J ^^ f {x)6{x — Xi) /\g' {xi)\, where the sum 
is over the simple roots of g{x), resulting in 

Taking — cq ^ C3 < cq as the integration boundaries yields the estimate [62, 63] 

^(A^e)~-^log(co), (2.14) 

where cq is small. The crucial feature is 

P(7Ve) oc 7V-3 , (2.15) 

indicating that shorter bursts of inflation are exceedingly more likely than longer ones ^ . This 
general feature has been observed in the two dimensional string landscape of [63] as well as 
in [59, 60, 79] and Sec. 4, using the toy landscape of Sec. 3. We would like to emphasise that 
both, the derivation and the result above, are not new and can be found in [62, 63]. 

Such a suppression raises an important question: why does our universe require around 
60 e-folds of inflation instead of say 40 (which would be more likely) to be consistent with 
observations? Based on the requirement that gravitationally bound structures form, an 
anthropic bound on (negative) curvature, and thus the minimum amount of inflation needed, 
was derived in [62], yielding N structure > 59.5 (see also [81, 82] for earlier work), close to the 
observed value. 

3 A Softly Bounded Potential 

From the brief review in Sec. 2 we learned that saddle-points are more likely than minima 
or maxima in a random potential (in the absence of upper/lower boundaries), the more so 



®If inflation does not take place near a saddle, but a shallow linear slope, the scaling becomes P{Ng) oc A^^* 
[62]. It should be noted that this result depends on the parametrization of the linear slope in [62] and the 
choice of a flat measure for those parameters. This freedom of choice in the measure is unavoidable in the 
phenomenological bottom up approach we employ in this study. 



the higher the diniensionahty of field space D. Further, for D 3> 1, inflation is most hkely to 
take place near a single saddle-point (if at all) and last for the shortest duration compatible 
with observational constraints. 

Here, we would like to investigate how bounds on the potential affect the ratio of saddle- 
points to extrema. Let us start with the simple potential proposed in [58-60, 79] 

V 
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where Xi = ipi/Mp, Mp = (SvrG) ^'"^ = 1 and aj^„jj^,bj^,„jj^ are Gaussian random variables 
with zero mean and variance 

^^i...^.=«-p(-i:^)' (3-4) 

and we defined 

D D 

Vji,...,Jd = «Ji...Jd cos^ JjXj + feji...jo sin^ JjXj . (3.5) 

j=l i=l 

In essence, the potential is just a truncated Fourier series with suppressed coefficients 
for higher wavenumbers Jj ^ 1, such that inflation can be accommodated [59, 60, 79]. 
Ji = • • • = Jd = is omitted to keep the potential levelled around V^ = 0. Opposite to 
a full Fourier series, only positive integer values for Jj are considered in (3.2) to keep the 
number of terms manageable as D increases; this comes at the price of sacrificing spherical 
symmetry, but will have no qualitative effect onto the results of this paper . Further, to 
ease comparison with [59, 60], we also set all coefficients to zero if a single index vanishes, 
Jj = (this choice is not physically motivated, but does not affect our conclusions either). 
Below the periodicity scale but above Ax ~ l/\/Dn ^, such a potential is certainly random, 
reminiscent of a beginners sky-slope, i.e. without cliffs or sharp bumps. To ease notation, 
we use 

J = Ji,.. . ,Jd , (3.6) 

X = xi, . . . ,xd , etc., (3-7) 

in the following. 

If points Xa and xj in field space are separated by more than Ax = |X(j — x^l ~ l/V-Dn, 
the value of the potential at such distinct points is again a Gaussian random variable with 

^AU analytic estimates in this paper could easily be done with the full Fourier series, but numerical 
comparisons would become much more challenging. 

^Coefficients of modes with wavelength below l/VDn are strongly suppressed, see (3.4); thus, the potential 
is not random on small scales. 



zero mean. The easiest way to see this is to consider first the same point Xc in different 
reahzations of the random potential (i.e. consider an ensemble): the Vj are the sum of 
two random, Gaussian variables, aj and bj, with zero mean, variance according to (3.4) and 
constant coefficients (the cosine and sine factors). As a consequence, the Vj are also Gaussian 
random variables with zero mean and the same variance 



av. 




dVjV fdVj^^ 



where the location Xc dropped out since 



+ (^) 1=^J' (3-8) 



dVjV fdVj^^ 




dajj +V^' =(cos>^J.xn +(sin>^J,xn =1- (3.9) 

Therefore, the value of the potential at a given point 

y(x,) = ^Fj(Xe) (3.10) 

J 

is also a Gaussian random variable with zero mean and variance 

/ DmT\ ' 



<yl 



E-i-T ^ (^"' 



To evaluate the sum, we first plugged in (3.4), approximated the sum by an integral which 
we could then evaluated in the limits (n + \)^^2/{Dn) » 1 and yj2/{Dn) <^ 1, that is for 

n>i:>. (3.12) 

Let us now turn our attention to distinct locations in the same random potential: as long 
as points are separated more than l/^/Du (but less than the periodicity scale), the corre- 
sponding values of the potential become uncorrelated; however, since the coefficients aj and 
6j are chosen only once, the potential is strictly bounded from above and below (sines and 
cosines are bounded); nevertheless, for \V\ <^ Sj o^J ~ {Dmr/4)^''^ we can approximate 
the probability distribution of F by a Gaussian in the limit of large n and D (invoking the 
central limit theorem); we may further approximate the variance of this Gaussian by (3.11) 
obtained from different realizations of the potential at the same pointy that is obtained by 
considering an ensemble of potentials^. 

This reasoning based on an ensemble is identical to the one used in [59, 60]: to in- 
vestigate inflation, the same starting point x = was used in different realizations of the 
random potential, instead of sampling different initial conditions in a single realization of the 
potential. The use of ensembles is of course nothing new and lies, for instance, at the heart 
of statistical mechanics. We shall use an ensemble of potentials again in Sec. 3.1 to compute 
properties of the Hessian at critical points. 

In the limit of large Dn, l/^/Dn approaches zero, indicating that the value of the 
potential at any point is an independent random variable. Obviously, such a potential would 



®Note that J]]j aj /av ~ {Dmr/2)^'^ 3> 1, providing a wide validity range. No important results in this 
paper depend on the actual value of ay in (3.11), and we treat ay as a free parameter in later sections. 



not allow slow roll inflation to take place and is therefore of little interest to us. Thus, we use 
n = 5 for all numerical computations in this paper, while treating the value of the potential 
at distinct critical points as independent, Gaussian random variables with zero mean and 
variance according to (3.11). By distinct we mean that critical points are separated more 
than 1/V-Dn; if they happen to be closer, we are essentially dealing with the same flat region 
in field space, rendering a distinction physically meaningless. 

An important property of our potential is that it is well approximated by a softly 
bounded one: we expect saddle-points to be ubiquitous for \V\ <^ ay-, in line with the 
conclusions of Sec. 2.1, whereas minima/maxima should become more likely for \V\ ^> ay^ 
which we shall prove now. 

3.1 The Hessian 

We denote the value of the potential in (3.2) at a critical point x by V{yi) = V, that is, 

n 

j=i, 

which, according to (3.11), is a Gaussian variable with variance ay ~ {Dmr/S) ' . We need 
to enforce this relationship if we want to retain the boundedness of the potential; to this 
end, we chose to replace l^i,...,i(x) by V and the shortened sum above. Hence, the Hessian 
becomes 

n 

= -V-Y,i'^kJi-l)Vj (3.15) 

j=i 
= -{V + hki), (3.16) 

for k,l = 1, . . . ,D; the sum covers the full range again, since Ji = • • • = Jd = 1 does not 
contribute. The hki are Gaussian distributed with zero mean and variance 

(3.17) 



'kl 



a^/3 for k^l, 



where we replaced sums by integrals, which we evaluated again in the limit n^ D. It should 
be noted that even though the h^i are not independent, we will treat them as such in the 
following, based on universality in random matrix theory. ^'^ 

Thus, given a critical point x with V'(x) = V, the Hessian is a symmetric matrix of 
Gaussian distributed real numbers with mean —V, variance a on the diagonal and cr/^/S 
otherwise. The novel feature is the presence of a non-zero mean for matrix entries which 
encodes directly the presence of the soft upper and lower boundaries in the potential, which 
entered by enforcing the constraint V'(x) = V. 



^"Correlations between the hki arise due to the limited number (2(n^ — D)) of terms in the potential, as 
well as the constraints VV|x=x = 0, which we do not impose. This harsh approximation becomes good if the 
potential contains a large number of independent Gaussian parameters, that is in the limit n^ D. 
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Figure 1: The probability Pmin of a positive definite Hessian (a critical point is a minimum) 
is computed for an ensemble of random matrices with Hij = —V — hij for {a) D = 2 and 10^ 
matrices per bin, (b) D = 2> and 5 • lO'^ matrices per bin, (c) D = 4 and 10^ matrices per bin, 
(d) D = b and 3 • 10^ matrices per bin; hij = hji are Gaussian random variables with zero 
mean, variances an = a on the diagonal and aij = cr/^/S for i ^ j, to mimic the Hessian 
of our toy random potential in (3.2). The dotted line is a two-parameter fit based on the 
Ansatz in (3.20), resulting in 0,^^ and p^un^ in Tab. 1, where we also compare Pmin{—oo) to 
the upper bound in (3.22) with (3.21) and p = 0.6. The D = 5 case appears noisy due to 
Pmin ~ 10~^ while drawing only 3 • 10^ matrices per bin, see Fig. 2. 



Since we are primarily interested in the effect of V, let us consider an even simpler 
Hessian by ignoring the off diagonal elements when computing eigenvalues, which become 

Xk = -{V + hkk). (3.18) 

Consequently, the probability that a critical point at height F is a minimum becomes 



Pmin =U{l- erf {v/iV2a)')U ; 



(3.19) 
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Figure 2: The probability Pmin as in Fig. 1 (d) for D = 5 and smaller values of V (a new 
drawing of random matrices is used); we observe that any features below V/a < —4 are 
indeed only due to noise (as expected), caused by drawing only 10^ matrices per bin. 

the probability for a maximum follows by replacing V by —V, and the chance for a saddle- 
point is Psad = 1 — Pmin — Pmax- For F = 0, we recover the result of (2.3), P = (1/2)^, 
that is saddle-points are the most likely occurrence. However, for large negative values of the 
potential, V < cr\/2erf"^(l - 2^~^/^) ss -few x a (i.e. for D ~ 10^), the probability that a 
critical point is a minimum approaches one. 

For the particular potential in (3.2) the variance scales as cr oc D '^Dn, so the range 
of ubiquitous saddle-points increases rapidly with increasing D, but if a different scaling of 
the coefficients' variance crj were chosen, one could alter this scaling easily. To keep our 
conclusions in subsequent sections as general as possible, we treat a as an unspecified free 
parameter from here on. 

Before we discuss consequences for inflation, lets step back and consider the more gen- 
eral, non-diagonal Hessian in (3.16); in this case, one could go the arduous route of computing 
the joint probability distribution (JPD) of eigenvalues and investigate the probability of a 
large fluctuation leading to a positive definite Hessian, as in [45]. However, a simple Ansatz 
(tested numerically) suffices for our purposes: as we saw in Sec. 2.1, the probability at y = 
is quite universally of the form P oc exp(— cD^) with c = 0(1) and 1 < p < 2 for JPD not 
too different from the Gaussian unitary ensemble. In the limit \V\ <C o", we could use the 
result of [46] in (2.5) if y <C 0.1, to investigate at which D the Gaussian suppression starts 
to dominate, i.e. when Dc2/ci > 1; however, for our particular potential the ratio of the 
variances is y = l/\^, naively indicating that the Gaussian suppression oc e~^ dominates 
for any value of D, but the analytic estimate in (2.5) is of course not applicable any more. 

To heuristically add the effect of a softly bounded potential, we make the Ansatz 



''iWa 



DP 



(3.20) 



with p and c as free parameters (specific to the Hessian at hand). For V <^ —a, the probability 
is bounded from above by the limiting value -P^™^* ~ {2e~'^)^'' . For D = 1,2 one can easily 
compute this limit analytically, leading to C£)=i = ln(2) and C£)=2 = lii(2)(l + 1/2^). We 
tested (3.20) numerically for D = 2,..., 5, see Fig. 1 and Tab. 1, finding good agreement 
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Table 1: Comparison of the Ansatz in (3.20) with Cu^n in (3.21) and p 

upper bound on P^^ in (3.22)) to an ensemble of random matrices with Hij — —v — ,oij, 



0.6 (leading to the 
-V - hi 
variances an = o" on the diagonal 
jj — cr/\/3 for i y^ j. See also Fig. 1. As expected, the probability of encountering a 
minimum is even more suppressed than Pi™'* , since the Hessian is far from diagonal. 



hij 
and a. 



hji are Gaussian random variables with zero mean, 



Dimension 


'-'nuin 


Cjimit 




/'nuni 


^™n(-00) 


plimit/ A 


1 


- 


ln(2) « 0.69 




- 


- 


1 


2 


1.17 


1.15 




0.61 


0.48 


0.5 


3 


1.69 


1.41 




0.59 


0.14 


0.25 


4 


2.20 


1.60 




0.60 


0.033 


0.125 


5 


3.12 


1.75 




0.52 


0.0038 


0.0625 


D:$>1 


- 


{1 + D^-P-D' 


'P) ln(2) 


- 


- 


2-D+l 



between the shape of our Ansatz and the full numerical result if we choose p 
by setting the normalization factor to^^ 

D-1 



c = ln(2) 1 + 



DP 



0.6. Further, 



(3.21) 



we arrive at an upper bound of encountering a minimum 



'(V 



-oo 



,--D+l 



(3.22) 



which is compared with the numerical values P^i^i^ ~^ — oo) in Tab. 1. One could improve 
the fit by including two exponents pi and p2, taking over some of the D-dependence of c 
at the expense of a second parameter, but we would like to keep the Ansatz as simple as 
possible (with the above c the exponential suppression is underestimated). Higher values of 
D become hard to treat numerically, due to the strong suppression of Pmin', hence, we use 
(3.20) with the above c for analytic estimates in the next sections, particularly to derive a 
lower limit on the magnitude of the expected contribution to the cosmological constant after 
an inflationary trajectory settled in a local minimum. One should keep in mind that the 
parameters are model dependent, but the (super)exponential suppression of Pmin{—oo) is 
generic. 

3.2 Where are the Minima? 

Since minima become exceedingly rare as D is increased, one may wonder at which height 
of the potential most of them lie. We will first give an upper bound of this height based on 
the results of Sec. 3.1, which we subsequently compare with numerics. 

Consider our toy potential in (3.2), which was chosen to be relatively smooth by sup- 
pressing higher frequencies with the aim to aid inflation. What is the number of critical 
points ric that are encountered if traversing a potential difference of order AV ~ a? Noting 
(I W|) ~ (T, we deduce nc ~ 1. We keep ric general for the time being, keeping in mind that 
it is a characteristic of the potential ^^. 

^^This is motivated by analytically investigating eigenvalues for D = 3 in the limit V <S^ —a. 
^^We consider classes of potentials for which ric is independent of D. 
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Figure 3: The probability of finding a minimum at height V is plotted over V/a for the 
potential in (3.2) and D = 1,2,3: we generate 10000 realizations for each D and search for 
minima according to the algorithm in Sec. 3.2. The dotted line is a Gaussian fit, whose mean 
is compared to the analytic upper bound of (3.26) in Tab. 2. 

Since the error-function is a bit cumbersome to deal with, we approximate Pmin in 
(3.20) by a step function 



where 



PrmniV) « 2^-^e(y, " V) 



y, = aV^erri(l-e^(2-^/^^)) 



(3.23) 



(3.24) 



and we used c from (3.21); Vc is chosen such that Pmin(yc) = -fmTn (~°o)/2. The probability 
of not finding a minimum when moving down the potential can thus be approximated by 



nn{V) 



[ 1 irr 



-oo))- 



-iV-Vc)nc/a- 



(3.25) 



for y < T4; if the distribution of the number of minima versus height is close to a Gaussian, 
which is the case here, we can estimate the position of the peak Vp by setting Pnominiyp) = 



14 



Table 2: The peak position of the histograms in Fig. 3 is compared to the analytic upper 
bound in (3.26) with ric = 1 and for different dimensions of the field space D. The histograms 
are the result of searching for minima in 10000 realizations of the potential in (3.2). 



Dimension 



v;"""Vcj 



T^limit 



1^ 



1 


-1.1 


-1 


2 


-2.2 


-1.5 


3 


-5.1 


-3.1 


Z)>1 


- 


-2^-Mn(2) 



1/2 and solving for 



K = v; + 



(J 



ln(2) 



ncln(l- 2-^+1) 



(3.26) 



For large D we get V^ ~ —a1 ln(2)/nc as the leading order contribution. For D = 10, 
this approximation is already better than one percent. Since we used the upper bound for 
the probability of finding a minimum in (3.20) with (3.21) and p « 0.6, the above provides a 
conservative upper limit on the height where most minima are encountered and thus a lower 
bound on the magnitude of the contribution to the cosmological constant. 

To check this bound, we generated 10000 realizations of the potential in (3.2) for D = 
1,2,3 and n = 5. To mimic the search of a minimum during inflation, we evolve the field 
equations in an expanding background^'^ , see Sec. 4, starting from x = 0; we lift the potential 
up by 5cr to keep Hubble friction high. Whenever the field comes to rest, we note the height 
of the un-shifted potential, and repeat. 

The analytic bound of the peak position is compared to the one in the histograms 
(Fig. 3) in Tab. 2; we observe that Vp does indeed provide the desired bound according to 
the limited range of D that are directly accessible for comparison. We further observe that 
the variance of the histograms is comparable to the location of the peak, 



0"n 



o{%. 



(3.27) 



which is expected given the gentle increase of the full Pmin in (3.20) from zero to Vp. 

As a direct consequence of the exponential suppression of Pmm(— oo) = 2^^+^, the 
position of the peak V^ shifts rapidly to increasingly negative values as D increases, even 
for our upper bound: for as little as ten fields we already have t^(10)/(T k, —356 and for 
field numbers commonly encountered in the string landscape, D ~ 500, we find t^(500)/(T ~ 
_2q150^ Thus, if inflation were to take place in such a potential, the ultimate resting place 
of the fields would most likely be in contradiction to a positive cosmological constant, even 
for extremely small a: 

Opposite to our naive expectation, the likelihood of finding a positive cosmological constant 
decreases rapidly as D increases - the huge growth of the absolute number of critical points 

is irrelevant. 



This method for finding minima is biased towards lower values of the potential; however, this search 
algorithm yields exactly those minima relevant for the vacuum selection during infiation. 
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This is one of the main results of this paper which we test numerically in Sec. 4 for low D. We 
comment further on the consequences for vacuum selection and the feasibility of multi-field 
inflation in general in Sec. 5. 

We would like to stress again that we severely overestimated the probability of finding a 
minimum; as a consequence, the above problem of encountering a large negative contribution 
to the cosmological constant is more severe in a realistic landscape. 

4 Inflation in a Random Potential 

We generate ensembles of random potentials defined in (3.2) with n = 5 and Z) = 1,2,3 to 
investigate the feasibility of inflation and its expected duration. The equations of motion are 
the Friedmann equations in a homogeneous and isotropic universe (A; = 0), 

3iz^ = y + X:f, (4.1) 

and the Klein-Gordon equations in an expanding universe 

dV 

ipi + 3H^, = -^—. (4.2) 

We are interest in the amount of inflation as measured by the number of e-foldings 

iV= f H{t)dt. (4.3) 



4.1 A Cut-off for large D 

Since the number of terms in the potential increases fast even for our already truncated 
Fourier series in (3.2), 

A^P,, = 2(n^ - £)) , (4.4) 

performing many runs becomes computationally challenging. To ameliorate this technical 
problem, one can truncate the series by setting coefficient oj and 6j to zero if their magnitude, 
chosen randomly according to (3.3) with variance (3.4), is below a prescribed cut-off. As a 
consequence, the potential becomes smother; hence, one should be careful to choose the 
cut-off not too high in order to preserve the qualitative features of inflationary regimes. 
Practically, we would like to keep a small fraction of coefficients only, F < 1, which we can 
identify on geometric grounds with 

Y jD 

IP _ ^D'Jmax /^ c\ 



2(n^ - D) 
where 



J= Ul 



\ 



D 



E J" ' (4-6) 



Vd = vr^' ^/r(L'/2 + 1) is the volume of the unit sphere in D dimensions, F is the Gamma 
function (r(x) ~ y/2TTx^~^''^e~^ if it is expanded for large arguments) and we assumed 
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Jmax < n ^'^- Since the variance in (3.4) is a nionotonically increasing function of J, we 
expect that most coefficients with J > Jmax wih have values smaUer than 

CTmin = 0-J^a, (4.7) 

(4.8) 
2/D\ 

TiD/2 + 1)] . (4.9) 




How big can one choose (Tmin without drastically altering the evolution of fields? The 
energy scale in our potential is set by the coefficients with the lowest J and thus of order 
y ~ 1. If we want to discuss inflation compatible with observations, we would need to leave 
the dynamics untouched even if the kinetic energy K = ^^ ^i'^/2 is only a small fraction of 
V, i.e. slow roll suppressed K ~ eV ~ 10~^. Coefficients aj, 5j < -fC ~ 10~^ will at most lead 
to a hill/valley that is easily traversed, even on an inflationary slow roll trajectory. Thus, 
imposing a cut-off of order 

CTrmn < W'^ (4.10) 

should leave conclusions with regard to inflationary dynamics unchanged. 

We tested numerically the effect of amin = 0,0.04,0.1,0.15 for D = 1,2,3 onto the 
relative probability of flnding N e-folds of inflation, see Fig. 4, and found no statistically 
signiflcant effect for cut-offs as large as amin = 0.15. An explanation might be that more 
expansion requires a flatter potential, that is the absence of high frequency oscillations. 
Hence, all runs with sizable amount of inflation operate on potentials with aj and 6j close to 
zero for larger J; these coefficients are thus unaffected by increasing the cut-off to relatively 
large values. 

4.2 The Feasibility of Inflation 

Based on inflation near a saddle-point, we saw in Sec. 2.3 that the probability of inflation 
should decrease sharply with the number of e-folds, P oc A^^^, independent of the dimen- 
sionality D of field-space. This behaviour can indeed be observed qualitatively in our toy 
potential (3.2); the algorithm we use for counting the number of e-folds for a given D is 

1. Choose the coefficients aj,bj according to (3.3) with variance (3.4). If D >2, impose 
a cut-off amin = 0.14 and set all aj, 6j to zero if their magnitude is below amin- 

2. Choose the initial field values and velocities'^ as ipi = (pi = for i = 1, . . . , D. If 
V{0) < reject the potential and go back to step 1. 

3. Evolve the equations of motion (4.1) and (4.2) numerically and compute the number 
of e-folds in (4.3). 

4. Stop the evolution if 1/ < or the fields get stuck in a minimum (i.e. no discernible 
movement within one e-fold). Save the number of e-folds up to this time and go back 
to step 1. 



^*For Jmax > n, the sphere does not fit into the cube with diameter 2n and we would count coefhcients 
already set to zero in the potential (3.2). 

^^Choosing the same starting point is justified since each run a new potential is generated. 
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Figure 4: The probability of A^-efolds of inflation in random potentials (3.2) with n = 5 is 
computed with and without imposing a cut-off amin = 0,0.04,0.1,0.15 onto the coefficients 
a J and 6j for D = 1,2,3 (top to bottom), see Sec. 4.1. We find no statistical significant effect 
for amin ^ 0.15, justifying the use of a cut-off. 

We plot the resulting probability for a large number of runs in Fig. 5 together with the 
best fit P = aN~^ yielding the values of a and /3 in Tab. 3. We observe somewhat smaller 
values for /3 than expected, which is not too surprising given that the analytical prediction 
ignores multi-field effects and assumes a head on encounter with a saddle-point, both of which 
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Table 3: Following the evolution in many potentials of the type 3.2 numerically and com- 
puting N according to the algorithm in Sec. 4.2, we plot the probability of finding N e- folds 
of expansion in Fig. 5. A power-law fit to P{N) = aN~" is performed; numerical values of /3 
are somewhat small than the theoretical expectation /3 = 3 in Sec. 2.3. a shows the expected 
strong suppression as D is increased. 



Dimension 


^ of runs cut-off amin 


a P 


1 
2 
3 


20000 

20000 0.14 

20000 0.14 


0.8 2.6 
0.3 2.3 
0.1 2.9 



become good approximations for large A^ only^^. This is in agreement with the conclusions 
of Frazer and Liddle [60] . In addition, they observed an increase of /3 if the potential becomes 
less featured, i.e. by increasing the horizontal mass ruh = ^i/xi (we kept nih = Mp fixed 
throughout). 

The values of a show a suppression with increasing D, which is expected, since it 
becomes less likely to encounter a flat stretch in field-space as D is increased. This feature 
was also observed in [60] for D = 1, . . . , 6. 

A novel feature not discussed in [59, 60] is the presence of a small bump in Fig. 5 at small 
A'^ ~ few. This feature is caused by a characteristic scale in the potential set by the modes 
with the lowest J; their effect can be understood qualitatively in one dimension: consider a 
single cosine as the potential, with unit amplitude and offset, V = 1 + cos(rE) (i.e. the J = 1 
mode only, with all other modes lumped into a constant); a randomly chosen initial value 
will neither be exactly at ^ = nor at y = 2. Thus, the probability of finding exactly zero 
or infinite amount of expansion is zero. Since P is continuous, there must be a maximum 
somewhere, which, since there is no other characteristic scale in the problem, should be near 
A^ ~ 0{1). The presence of more modes diminishes this effect, as evident in Fig. 5; however, 
since we strongly suppressed high frequency modes in the potential by choosing the variance 
according to (3.4), the bump prevails as D is increased from 1 to 3. Thus, we expect this 
feature to be characteristic for our type of "random" potential, but not a general prediction 
for random potentials. 

For D = 3, only 25 in 10000 runs ended up in a minimum above ^ = compared to 324 
and 1504 for D = 2 and D = 1 respectively (we only count runs that start with V > 0). Even 
though this is expected based on our analytic prediction that most minima have negative 
values of the potential as D increases, t^ < in (3.26) (see (3.20) with (3.21) as well as 
Tab. 1 and Tab. 2), we would like to stress the importance of this result, since a negative 
cosmological constant is in contradiction with observations. 

Observables in runs with more than fifty-five e-folds of inflation, i.e. in runs compatible 
with our universe, would be straightforward to compute; however, since almost all runs 
end up in an ADS vacuum (incompatible with observations), and our toy potential was 
not theoretically motivated anyhow, we refrain from following this avenue. If one ignores 
the negative cosmological constant, one may generate enough runs with sufficient inflation 

^^Furthermore, the smaU bump at low A'^ (see discussion below and Fig. 5) is not well fitted by a power-law, 
leading to a decrease of /? in a simple regression over all runs; this problem could be avoided by performing 
substantially more runs and including only the large A'' tail beyond this feature. 
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Figure 5: The probability of A^-efolds of inflation in random potentials (3.2) with n = 5 
is compared to a power-law P{N) = aN~^ with coefficients in Tab. 2.3 for D = 1,2,3, 
according to the algorithm in Sec. 4.2. For each D we perform 20000 runs. The relative 
importance of the bump like feature diminishes as D is increased. 
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for D ~ few (with modest computational effort) to provide statistics on observables such 
as the scalar spectral index, the tensor to scalar ratio or Non-Gausianities (NG); this was 
done in [59, 60] for D < 6: as expected, since almost all acceptable runs include inflation 
of the slow roll type, Ug and r are usually compatible with observations (in the sense that 
|ns — 1| <C 1 and small, but potentially observable r). Frazer and Liddle observe a dependence 
on the horizontal mass mh (less spread for increased rra^), but essentially no dependence on 
the dimensionality of field-space D. This is again expected, since the inflationary regime, 
taking place near a saddle-point, can usually be described by an effective- fleld-theory of 
much lower dimensionality. Consequently, NG should also be small, with two caveats: flrst, 
NG can be amplifled temporarily by a sharp turn in field-space (such turns are rare for 
the smooth potentials under consideration); second, as D increases, it becomes less likely 
that an adiabatic regime had been entered at the time of horizon crossing, rendering multi- 
field modes essentially unpredictable with regard to NG (they continue to evolve [41, 42]"*^^). 
Nevertheless, the expectation of commonly suppressed NG based on the horizon crossing 
approximation is consistent with the full numerical results in [60]. 

To conclude, even though dynamics in random multi-field potentials can be complicated, 
the overall features, such as the likelihood of inflation or the vacuum state after inflation, 
can be inferred from simple arguments in random matrix theory: as D increases, inflation 
becomes increasingly unlikely, is of short duration if it occurs, and ends up most likely in an 
ADS vacuum if the potential is bounded and symmetric around zero. Consequences of the 
last observation are wide reaching, as we shall see in the next section. Further, if inflation 
takes place, the inflationary regime can usually be described by an effective-field-theory 
of low dimensionality, with all of the observational consequences this feature entails (non- 
Gaussianities are an exception, since a prediction for i.e. CMBR fiuctuations requires the 
knowledge of the evolution of the entire field content until the universe entered an adiabatic 
regime) . 

5 Comments on the Vacuum State after Inflation 

5.1 Is Multi-Field Inflation unlikely? 

After inflation near a saddle-point in a higher dimensional landscape, the flelds evolve rapidly 
since the encounter with another saddle is unlikely [61]. We saw in Sec. 3 at the example of 
a particular softly bounded potential that minima are more common once V ^ Vp <^ —a, 
see eq. (3.26), that is, once the height of the potential is at a rare low value. Dynamics drive 
flelds to such rare locations, as we saw in Sec. 3.2 and 4. We argued that this feature is 
generic for all softly bounded random potentials. The distribution of minima, and thus for 
values of the cosmological constant, is far from even around zero, opposite to the assumption 
usually made to address the cosmological constant problem by means of anthropics coupled 
with eternal inflation, see i.e. the review [113]. If the value of the (metastable) vacuum after 
inflation is to be positive, as in our universe, we are forced to conclude one of the following: 

• Directly after inflation, the effective dimensionality of fleld-space is low so that an 
exponential suppression of the probability to encounter a minimum is absent; as a 
consequence, a landscape with an equal likelihood of positive and negative potential 
values can still be consistent with a small but positive cosmological constant. This 



^'^At worst, one needs to follow their evolution up until the universe is reheated. Unfortunately we lack a 
proper understanding of (p)reheating in most multi-field models. 
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requires the decoupling of almost all (massive) moduli, which may be reasonable as 
argued in [46]. 

• If the effective dimensionality of field-space remains large after inflation, the most likely 
vacuum to be reached dynamically is at rare, low values of the potential. For a positive 
cosmological constant to be likely, the potential on the landscape needs to be strongly 
yet precisely shifted towards positive values (or sharply bounded from below at F = 0). 

The second option requires our toy random potential in (3.2) to be shifted upward by many 
standard deviations, of order a2^^^ ln(2). Such a potential does not appear to be natural in 
string theory. We are thus left with the first option, which in turn raises another question: 

If the effective dimensionality of field- space directly after inflation is low, why should it be 

large during/before inflation? 

It appears highly fine tuned that after the (rare) encounter with a saddle-point (as needed 
for inflation in a higher dimensional random potential) the dimensionality of the landscape 
suddenly changes from D ^ 1 to D^ff ~ few. This is akin to running into a four star 
hotel in the Amazonian jungle (after being randomly dropped in said jungle), and after a 
relaxing holiday, one finds a convenient highway leading straight back to ones' porch (without 
encountering any junctions etc.). Without invoking anthropics (see Sec. 5.2), this needed 
tuning is a strong indication that only a few light fields are dynamically relevant during 
inflation^'^; hence, even though multi- field models appear to be natural in string theory with 
its plethora of potentially dynamical moduli fields, branes and fluxes, such scenarios are 
strongly disfavoured by the observation of a positive cosmological constant. A much more 
compelling scenario uses just one light field, so that an effective-field-theory can be applied 
[93-95], (see [83, 96-101] for selected applications)^^. 

Our conclusions are based on the purely classical evolution of scalar fields, invoking 
conservative upper bounds for the probability of encountering a minimum. The inclusion of 
tunnelling would strengthen our claims, since the added mobility [103] makes trapping in a 
metastable de-Sitter vacuum even less likely (the effect of resonant tunnelling [104] and giant 
leaps [105-107] improves mobility greatly as D is increased). 

However, our inability to reconcile the smallness of the cosmological constant with 
expectations from quantum-field theory is a strong indicator that we are still missing some 
crucial physics and we should be cautious to make such bold claims; hence, we do not 
advocate to throw out entirely multi- field models based on the above reasoning; nevertheless, 
to explain why we got stranded at such a peculiar vacuum after inflation is an additional 
challenge that needs to be addressed. 

It should be noted that the random potentials we focussed on show the same degree 
of "randomness" in all directions, opposite to landscapes leading to a meandering inflaton, 

^^The resulting presence of many, necessarily heavier moduli fields during inflation can still lead to inter- 
esting observational signatures, see i.e. [83-92]. 

^®By single-field EFT we mean that at any point on the inflationary trajectory only one direction has a 
light mass; the resulting trajectory can still be strongly curved (the corresponding potential would need to 
contain a gorge along which the inflatons are channelled). This notion of EFT is more general than the static 
definition of heavy and light fields used in [102], and still allows for the observational signatures in [83-92] 
that are absent in the treatment of [102]. 
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which arises if a small scale random potential is overlaid onto a smooth slow-roll slope, as in 
[27, 28]; for the latter our conclusions don't apply. 

Further, multi-field models with simple, similar potentials for all fields, such as in infla- 
tion driven by axions [16]), tachyons [17, 18], M5-branes [19-21] or D-branes [22-26] etc., are 
not affected by our conclusions either. However, these assisted inflation scenarios [108-111] 
are not the result of sampling all possible models in string theory, but a consequence of a 
deliberate choice by model builders. Thus, even though these models are concrete and pre- 
dictive, they tell us little about the models and predictions that are likely in string theory. 
More complicated setups, for which our conclusions are relevant, seem to arise naturally in 
string theory. 

Lastly, we would like to comment on the motivation to invoke several light fields during 
inflation: all current observations are fully consistent with the simplest models of single-field 
inflation. Concretely, recent progress on the inverse problem of reconstructing inflationary 
dynamics from the data by means of Baysian statistics [31-33] (see also [34-37] for general 
model discrimination), requires only two non-zero slow roll parameters [33], and a simple 
single-field potential suffices. If primordial non-Gaussianities were observed, we would be 
forced to include more complicated dynamics, such as a sudden turn or a temporary violation 
of slow roll among many other possibilities; however, most of these mechanisms do not require 
any additional light fields. 

5.2 A Way out via Anthropics? 

It appears conceivable to evade the conclusions of the last section by means of anthropics in 
an eternally inflating universe (see [112, 113] for reviews). How would such a resolution look 
like? 

First, one employs anthropic bounds for both the upper [55] and lower [56] values of 
the cosmological constant (CC). If a is much bigger than the range of anthropically allowed 
values of the CC, the probability distribution in (3.20) is essentially flat in this interval (in 
our concrete toy model the variance a in (3.11) increases sharply with D, in accord with this 
argument), even though it is not flat over the interval of all possible values; hence, ending 
up in a deSitter vacuum in this interval is roughly as (un)likely as finding oneself in an ADS 
one, which is encouraging. 

But how did we get there? Anthropic bounds on negative spatial curvature as well as 
observations of a nearly scale invariant scalar power-spectrum require around 60 e-folds of 
slow roll inflation [62], so our universe must have encountered a saddle-point or extremely 
shallow region in its past. The hight of this saddle needs to be large enough to allow for 
subsequent reheating and nucleosynthesis, many orders above the current value of the CC 
[113]. 

These events, finding such a saddle-point, having more than ~ 60 e-folds of inflation 
and ending up in a low vacuum with positive energy, are extremely unlikely, the more so 
the larger D is. Thus, a mechanism is needed to sample many trajectories, such as eternal 
inflation. Since minima are rare for large V, it appears reasonable that eternal inflation is 
driven by fields high up in the potential (kept there by quantum fluctuations); however, our 
universe could also arise due to fields being trapped in a metastable vacuum prior to slow roll 
inflation - after all, once started, eternal inflation populates all vacua [54]. The latter yields 
an open universe after tunnelling [114-116], leading to the already mentioned anthropic lower 
bound on N [62]. 
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Once eternal inflation is invoked, we run into a sever problem: all statistical arguments, 
such as the ones just made, become meaningless due to our inability to impose a mathemat- 
ically well defined measure on countable infinite sets [57] (see also [117]). A possible way 
out of this measure problem is to realize that eternal infiation is not past [118] but future 
eternal. Thus, at any given "time" the set of bubble universes in the multiverse is actually 
finite^*^. Then the conceptual problems discussed in [57] become more practical ones: how 
should we define a gauge invariant time in an inhomogeneous multiverse so that we only 
deal with finite sets, and at what time should we evaluate probabilities (taking the limit of 
time going to infinity is not appropriate)? The simplest choice, taking constant proper time 
slices, leads to the youngness paradox [58]. Other proposals, such as using scale factor time 
[119-121] (a global measure) or the causal diamond measure [122] (a local one), seem to work 
better (see [123] for a recent review of both global and local measures)^^. Particularly, the 
causal diamond measure apparently leads to the correct post-diction of the CC [125-127]. 
Nevertheless, how do we know which measure is the right one in our multiverse? In the 
end, it is our choice of slicing that determines observables, so in what way are we predicting 
anything? This unfortunate freedom of choice is present due to our (conceptual) inability to 
probe the multiverse. Maybe all we can hope for is to follow a pragmatic approach, choose 
a measure ad hoc such that some observables come out right and then make predictions for 
others. 

Following this line of thought, one needs to be aware that the entire inflationary dynam- 
ics appears to be set by anthropics. The resulting history of our universe, eternal inflation, 
the unlikely encounter with a reasonably hight saddle-point, inducing an unlikely long phase 
of slow roll inflation, followed by the extremely unlikely encounter with a metastable, low 
deSitter vacuum, is far from generic. 

As a consequence, studies that sample a random landscape, such as [59, 60], with the 
goal to extract distributions of observables (the scalar spectral index, the tensor to scalar 
ratio, non-Gaussianities etc.) are problematic if not all observational bounds are properly 
imposed. For example, in [59, 60] trajectories were considered regardless of the value of the 
CC. It is not clear to us that the presence of a metastable deSitter vacuum very close to the 
inflationary saddle-point has no bearing on observables; for example, non-Gaussianities are 
sensitive to dynamics until the universe reheated and an adiabatic regime is entered [41, 42]. 

One conclusion is clear to us: if the landscape in string theory is as complex as it 
appears to be (at the energy scales relevant for inflation), our universe's history is anything 
but generic. 

6 Conclusion 

We investigated the feasibility of softly bounded, random potentials in higher dimensional 
field-spaces for inflationary model building. Based on random matrix theory, we found saddle- 
points to be far more generic than minima or maxima for common values of the potential 
a D > few; as a consequence, inflation takes place most likely near a saddle, if it occurs 
at all. Further, shorter bursts of inflation are far more likely than longer ones, P oc N~^ 



■^°The common reason infinities arise in eternal inflation (and thus a measure is sought after), is the com- 
parison of all possible observers, future and past alike; however, such an approach violates causality: any 
predictions for our universe can only depend on events in the past, not possible events in the future. It 
appears that this fundamental principle of physics is contested in this field of research. 

^^We would like to point out that we disagree with some of the conclusions in this review, i.e. the "prediction" 
of the end of times [124]. 



24 



with /3 ~ 3. We tested these general expectations in a simple toy-landscape, chosen to be 
reasonably smooth while remaining softly bounded by means of a Gaussian factor {P oc 
exp(— y^/2(Ty)/(Ty), and we found good agreement. 

We followed with theoretical and numerical investigations of the final state after infla- 
tion in such a random potential. Since the probability of running into a minimum decreases 
dramatically for positive values of the potential as the dimensionality of field-space is in- 
creased, see eqn. (3.20), it becomes exceedingly unlikely to end up in a universe with a 
positive cosmological constant after inflation. We confirmed this expectation numerically for 
the same toy-landscape; for instance, in the presence of three fields less than 25 in 10000 runs 
(indiscriminate of whether they include inflation or not) ended up in a minimum above zero. 
Taken together with the observational requirement of around sixty e-folds of inflation, it is 
extremely unlikely that our universe originated from dynamics involving more than a cou- 
ple of light fields during inflation. To reach this conclusion, we used conservative estimates 
throughout and we deliberately avoided anthropic arguments, which could be used as a way 
out. 

As a consequence, the eff'ective-field-theory of single- field inflation in the presence of 
some heavier fields, see i.e. [83-86, 93, 94] for recent progress, could not only be entirely 
sufficient to tell us everything we will ever know about infiation, but appears to be favoured 
over more complex multi-field models in light of a positive cosmological constant. 
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